function ans=f1(t,T)
	global TOUT_M;
	global TD;
	global k;
	global ku;

	ans=k*(Tout(t)-TOUT_M)+(k+ku)*(TD-T);
end

